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Abstract 

In this paper, we continue our investigation of a one-dimensional, two- 
component, quantum many-body system in which like particles interact with 
a pair potential s{s + l)/sinh^(r), while unlike particles interact with a pair 
potential —s{s + l)/cosh^(r). For an equal number of particles of the two 
components, the ground state for s > corresponds to an antiferromag- 
net /insulator. Excitations consist of a gapless pair-hole-pair continuum, a 
two-particle continuum with gap and excitons with gap. For — 1 < s < 0, 
the system has two gapless excitations — a particle-hole continuum and a 
two spin-wave continuum. Using finite-size scaling methods of conformal field 
theory, we calculate the asymptotic expressions and critical exponents for 
correlation functions of these gapless excitations at zero temperature. The 
conformal structure is closely related to the Hubbard model with repulsive 

on-site interaction. 
64.60.Fr, 71.45Gm, 75.30Ds 
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I. INTRODUCTION 



We recently presented the exact solution to a one-dimensional, two-component, quantum 
many-body system of considerable complexity in the form of an asymptotic Bethe Ansatz 
calculation^ The two kinds of particles are distinguished by a quantum number cr = ±1, 
which may be thought of as either spin or charge. The system is defined by the Hamiltonian 

1 92 

H = - 9^+ E Vjk{xj-Xk), (1) 

where the pair potential is 

1 + ajak 1 — cToCrfe 



Vjk{x) = s{s + 1) 



2sinh^('x) 2cosh^('x) 



(2) 



We assume s > — 1. We call this the SC-model, for the sinh-cosh interaction. Thus for s > 0, 
like particles repel, while unlike particles attract. When like particles are near, the repulsive 
potential increases as while for large separations, both potentials decay exponentially 
with a decay length which we take as our length scale, and hence unity. The potentials 
might usefully be thought of as a screened potential. This system was first introduced 
by Calogeroi, who showed it to be integrable. Sutherlandil soon afterward showed that the 
system could be exactly solved, and gave the solution for a single component system. In the 
present paper, we continue our study of the SC-model with an investigation of the correlation 
functions and their critical exponents, at zero temperature, by methods of conformal field 
theory. 

The concept of conformal invariance in one-dimensional (ID) quantum systems at crit- 
icality constrains the possible asymptotic behavior of correlation functions and allows a 
classification into universality classes, distinguished by the value of the central charge c of 
the underlying Virasoro algebra.S For models with short-range interactions and a gapless 
excitation spectrum with a single Fermi velocity, we can determine both c and the critical 
exponents of correlators directly from finite-size corrections to the ground state energy and 
the low-lying exited states. In most cases, such models have been found to belong to the 
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universality class of the ID Luttinger model,! i.e. c = 1, and the critical exponents to vary 
as functions of the coupling constant of the corresponding conformal theory. 

Recently, various authors have extended these concepts to include multi-component sys- 
tem. with different excitation velocities, =ncl> a. the Hubbard model.l In general, one Bnd. 
a c = 1 Virasoro algebra for each critical degree of freedom, i.e. each gapless excitation with 
a unique velocity. It is then possible to construct the full theory as a semidirect product 
of these independent algebras. Again, critical exponents may be calculated from finite-size 
corrections but now they are functions of a matrix of coupling constants. 

In another recent development, the ideas of conformal field theory have been applied to 
models with long-range interactions such as the system.!'! It turns out that one can 
no longer simply read off the value of the central charge from the finite-size corrections to 
the ground state energy. However, one may still calculate the correct critical exponents 
of the asymptotics of the correlations functions from the finite-size scaling behavior of the 
low-lying excitations.0 

Our paper is organized as follows: In section ^ we briefly review the asymptotic Bethe 
Ansatz solution of the SC-model in the zero sector as obtained in Ref. |^. Section ^ outlines 
the Luttinger liquid approach for long-ranged models. We give arguments why the standard 
evaluation of conformal field theory for c fails for long-ranged models. For s > 0, there is only 
one gapless excitation corresponding to a single c = 1 conformal theory. For — 1 < s < 0, 
however, there are two gapless excitations with different velocities, so that we briefly review 



the main formulas for a semidirect product of two c = 1 Virasoro algebras. In section |V 
and section 0, we give expressions for the correlation functions and calculate their critical 
exponents from the Bethe Ansatz equations for the — 1 < s < and the s > cases, 
respectively. For simplicity, the two types of particles are assumed to be either both bosons 
or both fermions, although a mixed bose/fermi system can be studied along similar lines. 
We close our paper with section |VI] , where we briefly show how both the — 1 < s < and 
the s > cases fit together as s — 0''^. 
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II. THE BETHE-ANSATZ SOLUTION IN THE ZERO SECTOR 



Let us recall the results of Ref. 0: We have particles with a = —1 and N-^ with 
cr = +1, for a total of N = + N-^ and N-^ > N^. The zero sector corresponds to an equal 
number of up and down spins, i.e. A^^ = A^^^. For s > 0, pairs of up-down spins bind into a 
variety of bound states, or pairs, which we will label by m = 1, . . . , M{s), where M{s) is the 
smallest integer larger than s. Let there be of each type. Unbound particles correspond 
to spinons/ions in the spin/charge picture and there are Nq = N — 2J2i<m<M of these. 
Let us call particles with spin down spin waves; there are A^_i = Ni — J2i<m<M of these. 

Imposing periodic boundary conditions on the wave function and taking any particle, pair 
or spin wave around a ring of large circumference, yields the following set of Bethe-Ansatz 
equations. 

Lf]rnM-m, 27rl^(km) -|- ^ ^ ^ ^ ^m,m' Q^m ^m'\ 

-l<m'<Mk' , 

m = -1,0,1,..., M. (3) 

Here the I.m(k.m) denote the set of quantum numbers for each type of particle. Depending 
on the parities of A^^ and the particle statistics, the quantum numbers will be restricted 
to integer or half-odd integer values. Note, that for the spin waves, I_i ranges only over 
1, . . . , A^o- Om,m'(Mm — ^m') IS the phasc shift for the scattering of particles of type m with 
type m' and has been calculated previously.0 Note that 9mm'{k) = ~9mm'{—k) = 9m'm{k). 
Furthermore, we define 

0, m = —1, 

1, m = 0, (4) 

2, m = 1,2, ...,M(s) 

We can write the momentum and energy for a solution of as 

P= ri^^Ki, (5) 

-l<m<M km 

E = \ E v^HK.- E (6) 

^ -l<m<M k™ l<m<M 



n. 



m 



Here Km — s + 1 — m denotes the poles in the transmission and reflection coefficients that 
give rise to bound states of type m. 

For > s > — 1, there are no bound states and we will call this the unbound case in the 
sequel. We therefore have only two coupled equations for A^o particles with pseudo-momenta 
ko = {ki, . . . , /catq) and A^_i spin waves with rapidities k_i = (Ai, . . . , Ajv.J- 

+ E^JMk,-k), 
= 27rJ,(A„) + E?=!^-i,-i(Aa-A^) 
+ Efii^-i,o(Aa-A;,). 

The particle quantum numbers Ij and the spin-wave quantum numbers Ja are restricted by 
the parities of A^O) -^-i and the statistics of the particles to the following combination of 
integers and half-odd integers: If both spin-up and spin-down particles are bosons 

Ij = {No - l)/2 (mod 1), 

Ja = (A^-i - l)/2 (mod 1), (8) 

whereas for fermions, 

I- = N_i/2 (mod 1), 

Ja = (A^o + N^i - l)/2 (mod 1). (9) 

In the thermodynamic limit, i.e. L — > oo with fixed do = Nq/L, d-i = N^i/L, the 
ground state is a filled Fermi sea characterized by the distribution function p{k) of particles 
and (t(A) of down-spins. 

P(k) = ^ + ^I^cO'o,-^(k-^^)a(|,)d^, 

+ i^I^BOo,oik-h)p{h)dh, ^^^^ 
<7(A) = + ^/5c^-i,-i(A-/x)a(/.)d/. 

+ i^!^sOU,oi^-h)pih)dh. 
Here the prime denotes the first derivative. The values of B and C are fixed by the following 
equations: 




(T(A)dA = d 



p{k)dk = d[ 



-1 — 



- M 



(11) 



(12) 



where M. = {Ni^ — Ni)/2L is the magnetization per unit length. Let us now restrict our 
discussion to the zero sector, when A^^o = A^, A^-i = M. = Q and the hmit C of the 

spin wave distribution is oo. Then we can solve for the spin wave distribution by Fourier 
transform in terms of the particle distribution, which we then substitute into the particle 
equation, giving a single integral equation for the distribution of particles p{k): 



The excited states in the zero-sector are given by: (i) Remove a particle from the ground 
state distribution, and place it outside the limits; we call this creating a hole and a particle, 
and it gives a two parameter continuum, (ii) Remove a spin wave from the ground state 
distribution, and place it on the line with imaginary part equal to i; we call this creating 
two spin waves, one with spin up and the other with spin down. It gives a two parameter 
continuum of the type familiar from the Heisenberg-Ising model.0 Each of these two types 
of two-particle continua has a single Fermi velocity. Let us denote by vq the Fermi velocity 
of the first excitation and by f_i the Fermi velocity of the second. As has been pointed 
out in Ref. ^ the two velocities are in general not identical. The same is true of the 
Hubbard model with repulsive on-site interaction, and we will later make extensive use of 
the conformal results obtained for this model.iS 

For s > 0, which we call the bound case, the ground state in the zero sector consists 
of a spin fluid of type m = 1, and thus spin 0. In the ground state, the fc's for the pairs 
distribute themselves densely with a density t(A;), between limits ±D, normalized so that 




(13) 



Here the kernel 9' 



(k) is given as 




(14) 




(15) 
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The energy and momentum are given by 

P/L = 2 T{k)Mk = 0, (16) 

J-D 

E/L = r T{k)k'^dk - s^NjL. (17) 

J-D 

The integral equation which determines T{k) is 

I/tt = T{k) + ^ r 9'u{k - h)T{h)dh. (18) 

ZTT J-D 

The kernel of the equation, 0[i{k), is the derivative of the phase shift for pair-pair scattering. 

The low-energy excited states are given by the following: (i) Remove a pair from the 
ground state distribution, and place it outside the limits; we call this creating a pair-hole 
and a pair, and it gives a gapless two parameter continuum, (ii) Break a pair, to give two 
particles, one spin up and the other spin down; this also gives a two parameter continuum. 
However, there is a finite energy gap for breaking a pair. These are the spinous or ions, (iii) 
Excite a pair into a higher energy bound state, if allowed; these we call excitons, and they 
have simple single parameter dispersion relations. 

Let us denote the unique velocity of the excitations of type (i) by vi. The Bethe-Ansatz 
equations that describe these excitations may be written as 

N/2 

2Lkj = 27iHj{kj) + J2 - k)- (19) 

1=1 

Note here that kj is the pseudo-momentum of a pair, and is not the pseudo-momentum of 
an individual particle, which would be complex and of the form kj/2±is. The pair quantum 
numbers Hj are restricted by the parity of A^i, and bose and fermi statistics are given as 

Hj = (Ni + l)/2 (mod 1), (20) 

since pair-pair scattering is symmetric for pairs of bosons and pairs of fermions. The pairs 
will be singlets. 
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III. CONFORMAL APPROACH FOR CORRELATION FUNCTIONS 



A. finite-size scaling in conformal theories of Luttinger liquids 

The behavior of the correlation functions for a given one-dimensional model at large 
distances and low-temperatures is determined by the gapless excitations. Ill These gapless 
excitations are due to hydrodynamic fluctuations and it has been arguedi that the low energy 
physics of such a system may be described by the exactly solvable Luttinger model,lli the 
ID quantum version of the classical 2D Gaussian model. The Luttinger model is a critical 
system with continuously varying exponents and corresponds to the universality class of 
c = 1 conformal field theories.0 Application of conformal theory allows the calculation of 
these critical exponents purely from finite-size scaling arguments.lll 

The value of the central charge c may be read off from the following finite-size scaling 
formula 

TTV 

Eo ~ eoL - —c, (21) 

thus enabling an independent check of the above arguments. Here, Eq is the ground state 
energy of the finite system, eo is the ground state energy density in the thermodynamic limit 
and V is the Fermi velocity in the system. In short-ranged ID quantum models, including 
Bethe-Ansatz solvable models, the above universal picture is confirmed.^ However, for long- 
ranged models, straightforward application of this equation may lead to unphysical results.^ 
(we include the SC-model in this class, although its pair potential decays exponentially, 
since it can only be solved by means of the asymptotic Bethe Ansatz.) For instance, in the 
models c is predicted to be equal to the interaction strength, although independent 
calculations show that the critical exponents are those of the c = 1 universality class.0 
However, if one instead estimates c from the low temperature expansion of the free energy,0 
one does get the correct answer c = 1. 

We can, however, understand the failure of (B^) in long-ranged models. The crucial point 



is that, due to the long-range character of the interactions, finite systems will always "feel" 
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the particular boundary conditions chosen, so that ( pT]) includes an additional correction 
term Ehc, representing the boundary energy, and so 

Eo^eoL-—c+^. 22 

The low temperature expansion, however, uses boundary conditions instead for the time 
axis of the model and we thus have no such corrections. We may therefore write the free 
energy of a long-ranged c = 1 Luttinger liquid as 

F(T) ^ F(T = 0) - (23) 

Of 

Let us recall the main formulas for calculating the correlation functions and their critical 
exponents.0 Every primary field in a conformal field theory on an infinite strip of width 
L in the space direction gives rise to a tower of exited states. Let x = + denote the 
scaling dimension and a = A^ — A^ the spin of (f)±. Then the energies and momenta of 
these exited states scale as 

E{A^, N^) -Eor^—{x + N^ + N-), (24) 

P(A±, N^) -Por^'^{(T + N+ -N-) + 2Dkf. (25) 

Here A^"*" and are positive integers, v is the common Fermi velocity of the excitations 
and 2D is the momentum of the state in units of the Fermi momentum kf. Note that the 
quantities on the left hand side of these equations are evaluated with respect to the same 
boundary condition and therefore the above mentioned correction terms cancel. We may 
write the correlation functions of the primary fields at zero temperature (expressions for low 
but finite temperature may also be given) as 

exr)(—2iDkf) , . 

{*^.(x.t¥A.(0.0)> = ^— J-^^^iL_. (26) 

However, the excitation spectrum of the SC-model is quite different for the bound (s > 0) 
and the unbound (—1 < s < 0) case as we have argued in the previous section. Most 
importantly, the unbound case does not have a common velocity for all excitations anymore 
and so the formulas given above for a Lorentz-invariant conformal field theory can no longer 
hold. 
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B. conformal weights and the dressed charge in the bound case 



For the bound case in the zero sector, only the pair-pair-hole excitation branch is gapless. 
Thus there is only one excitation velocity and from the above arguments, we expect the 
dimensions of the primary operators to obey the formulae for a single c = 1 Gaussian 
model, i.e. 

A^iAN„D,) = ^(^D^^^±^J . (27) 

The coupling constant of this Gaussian model depends on the system parameters. It is 
sometimes called the dressed charge and may be calculated from the Bethe-Ansatz equations 
by means of an integral equation! 

ei(fc) = 2 + ^ r Uh)e[ Ah - k)dh, (28) 

where the constant is 2 because this excitation is a pair. However, we can also calculate 
^1 = ^i(-D) by purely thermodynamical arguments as follows: Let us change a given ground 
state configuration by adding pairs while keeping the Fermi sea at zero momentum, so that 
the excitation can be described by the pair (AA^^i, Di = 0). Then a second order expansion 
gives 

AE = -^i(AiVO + ^;^(AiVl)^ (29) 

where fxi = —-§§^ is the chemical potential for adding pairs and ki is the pair-compressibility. 
Comparison with (|2^) and (^7|) yields 

^\ = nviHidl = 7idi/vi. (30) 

In the last equation, we have used the well known relation vf = l/(/ti(ii). Therefore, by 
knowing the Fermi velocity of the pair - pair-hole excitations, we can calculate the scaling 
dimensions. 
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C. finite-size scaling and the dressed charge matrix in the unbound case 



For the unbound case, two excitation branches are gapless, giving rise to a particle-hole 
continuum and to a spin wave continuum, with Fermi velocities Vq and respectively. 
Thus, the finite size corrections of equations (^) and ( ]25|) now become 



27T 

E(AN,D) - ^0 ~ — 
L 



(31) 



P(AN, D) - Po ~ 1- AN^ Ti + N^ - + - NZ^ + 2DoA;/,t + 2(Do + 



(32) 



Here, the matrix V" = diag(fo, ti.i) and the excited state is characterized by the pairs 
AN = (AiVo, AA^_i) and D = (Do, /^-i). As before, A^^^ and A^^i are positive integers that 
label the descendant fields. The 2x2 matrix S is the generalization of the dressed charge 
and may be calculated by means of coupled integral equations. Thus if we denote the 
components of H by 



then 



eo,o(fc) 
^0,-1 (A) 
e-i,o(A;) 
e-i,-i(A) 



^ eo,-i(C) ^ 



1 + ll^B^o,o{h)0'o,oih-k)dh 

+ ^/5c^o,-i(/^)^-i,o(/^ - ^)d/i. 

+ ^J^^^o,oih)0'o,-i{h-\)dh 

+ 2^/5ceo,-i(/i)^^,-i(^-A)d/i, 

+ 2^/5ce-i,-i(/")^-i,o(/"-^)d/i, 

1 + lI^B^~i,o{h)e',^-i{h-X)dh 

+ lJ%i-,,-Me'_,^-l{^^-\)d^,. 



(33) 



(34) 
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Thus, the situation for — 1 < s < is analogous to the situation in the repulsive Hubbard 
model away from half-filling0'l and we may interpret equations (|3T| ) and ( p^ ) in terms of 
a semidirect product of two independent Virasoro algebras, both with c = 1. The scaling 
behavior of the energy and momentum in terms of the conformal weights and A^^ and 
the formulas for these weights as functions of the components of the dressed charge matrix 
S have been given in Ref. |^, and we will not repeat them here. The generalization of the 
correlation functions of the primary fields has also been given in Ref. 0. However, as before, 
thermodynamic arguments may be used to calculate the values of the dressed charge matrix. 

For the zero sector, i.e. = 0, the relevant equations simplify considerably. In this 
case, kf^i = kf^^ = kf = 7Tdo/2, and the dressed charge matrix S may again be expressed in 
terms of a single parameter = ^o{B), i.e. 

/ ^0 ^ 



(35) 



Thus the conformal weights A^ and are given as 



1 



± ^AiVo(2Do + D^i) + N^, (36) 

± ^(2AA^_i - ANo)D^i + N^^ (37) 

Note that the second equation is independent of ^o- However, there is an explicit dependence 
on the interaction strength s and only for s = do we recover the result of the Hubbard 
model. 

This s dependence can be understood by realizing that for the zero sector and — 1 < s < 
the Bethe-Ansatz equations of the rapidities k_i = (Ai, . . . , Xn^i) are essentially the Bethe- 
Ansatz equations of the Heisenberg-Ising model. The effect of the Bethe-Ansatz equations 
for the pseudo-momenta is simply a renormalization. Following Ref. |TT] we parametrize 
the anisotropy in the Heisenberg-Ising model by A = — cos(yu). Then the correspondence 
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is established by setting /i = —tts. Thus we may say that the behavior of the spin wave 
exitations changes from ferromagnetic at s ^ — 1^ (A = 1) to antiferromagnetic at s ^ 0~ 
(A = —1). Furthermore, we expect to see free spin waves at s = —1/2. This picture has 
been confirmed by a study of the transport properties of the SC-model which we present in 
another pubhcation.lli 

An integral equation can also be given for ^o, 

Uk) = ^ + ^ r Uh)e\h - k)dh, (38) 

where the kernel is as in equation (0). Alternatively, we may simply express in terms of 
thermodynamical response functions as 

^0 = TTVoKodl = ndo/vo- (39) 



D. correlation functions and conformal expansion 

Given the conformal weights, we now construct the asymptotic expressions for correlation 
functions. For — 1 < s < 0, we want to consider the following set of correlators: Let il!^{x,t) 
denote the field operator of a particle with spin a. Later, we will additionally restrict the 
statistics to be either bosonic or fermionic by restricting the possible values of the pair D. 
Then the field correlator — also called the one-particle reduced density matrix — is given 
by 

C4x,t) = {^^{x,t)^l{0,0)). (40) 

Defining the number operator n{x,t) = n-^{x,t) + ni{x,t), we write the density- density 
correlator 

Cn{x,t) = {n{x,t)n{0,0)). (41) 
The spin-spin correlation functions are 
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C:{x,t) = {S^{x,t)S\0,0)), 



(42) 



C^(x,t) = (5-(x,t)5+(0,0)), (43) 

where we used = {n^ —n^)/2 and = ip^ipi- Note that for systems that are rotationally 
invariant, such as the Hubbard model in zero magnetic field, these two spin-spin correlators 
are closely related, i.e. = 2C^. 

Following Ref. 0, we also consider the correlation function for singlet pairs, 

Cs^ng{x,t) = (^}(a;,t)^{(a;,t)^r(0,0)^i(0,0)). (44) 

Note that all these correlators are of the form {A{x,t)A^ (0,0)). By standard arguments of 
conformal field theory,@ we may deduce the leading terms and the critical exponents of the 
long-distance behavior of these correlators by expanding A in terms of the primary fields 
while minimizing with respect to D at the corresponding values of AN. Here the above 
mentioned restrictions on D will become crucial. This approach, however, will leave the 
expansion coefficients undetermined and at special points in the phase diagram, they may 
even vanish. 

For s > 0, the model exhibits a gap for breaking of pairs and there are no spin waves. 
Therefore the correlators (0), (0) and (|43|) will decay exponentially. Let us introduce the 
pair field operator \E'. The pair density - pair density correlator can be written in terms of 
the pair number operator p = as 

Cp{x,t) = {p{x,t)p{0,0)) (45) 

and the pair field correlator is given by 

C^{x,t) = {^^{x,t)'^ (0,0)). (46) 

As before, we can construct these correlators by an expansion in primary fields, minimizing 
with respect to ANi and Di. 
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IV. ASYMPTOTICS OF THE CORRELATION FUNCTIONS FOR THE 

UNBOUND CASE 



Due to the restrictions (P) and (P) on the quantum numbers of a given state, the numbers 
D = {Do, -D-i) are integers or half-odd integers depending on the parities of the pair AN = 
(AA'^o, AA^_i) and the statistics of ip^ip. In particular, for fermionic particles we have 

AiVo + AAU AiVo 
Do = (mod 1), D^i = (mod 1). (47) 

We can now apply the scheme for calculating the leading asymptotic behavior of the cor- 
relation function as outlined in the last section. Following our selection rules, we therefore 
have for a fermionic system 

: ANo = 1; AN_, = 1; 

Do = 0,±l,...; D_i = ±i,...; 
Cn : ANo = 0; AiV_i = 0; 

Do = 0,±l,...; D_i = 0,±l,...; 
AiVo = 0; AiV_i = 0; 

(48) 

Do = 0,±l,...; D_i = 0,±1,...; 
: ANo = 0; AiV_i = 1; 

Do = ±i...; I^-i = 0,±l,...; 
C,,„g : ANo = 2; AiV_i = 1; 

Do = ±|,...; D-i = 0,±l,... . 
This is identical to the results for the repulsive Hubbard model, and as in Ref. ^, we will write 
the critical exponents as functions of ^ = 2^q. However, there is an additional interaction 
strength dependence in the correlation functions due to the explicit appearence of s in 
equation (pTj). This is a novel feature and not true in the Hubbard model. It emphasizes 
the close correspondence of the Heisenberg-Ising model and the SC- model for — 1 < s < 
in the zero sector. 

Following the scheme outlined briefly in the last section, we calculate the leading asymp- 
totics of the fermionic field correlator in the SC-model to be 

15 



1 



Re 



X — ivot 



+ 



\x + ivot\^/'^+^'>/^^\x + zt;_it|H^V4(«+i) 
The density-density correlator is given by 

cos{2kfX + $1 



Re 



Aie 



.. fx + ivot 



X — ivot 



X — iv-it 



X + if _it 



X — iv-it 



+ A 



\x 



cos{AkfX + $2) 
\x + ivot\^ 



+A3 



x^ - (votf 



+ Aa 



X 



(50) 



[x^ + {voty]^ "[x2 + (7;_it)2]2' 

and since the selection rules for the density-density correlator are identical to the selection 
rules for the longitudinal spin-spin correlator, the above calculation holds for with dif- 
ferent constants and the replacement of Ai^ for ri^. Finally, for the transverse spin-spin and 
the single-particle correlator we find 

cos{2kfX + $) 



C^{x,t)^Ao- 



X + Wot\'^/'^\x + W_it\^^+'') 



1 



Csing{x, t) ~ Aq 



1 



Re 



Ai 



X + if _it 
X — iv^it 



1 



rRe 



X — ivot. 



(51) 



(52) 



Following equation (|39|) , we calculate ^0 from the Fermi velocity Vq. In Fig. ^ we plot the 
lines of constant ^0 in the {cIq, s) plane. Note that the value of ^(^0) at zero density is 
given by 2(1), whereas for finite densities and vanishing interaction strength s ^ 0^, we 
have ^ — s> 4 (^0 V^)- As expected, this is the same behavior as in the Hubbard model for 
vanishing on-site interaction strength u. In particular, the explicitely s-dependent exponents 
in the SC-model reduce to constant values as s ^ 0^ which are equal to the corresponding 
exponents in the Hubbard model. However, we can not bound 6 between those two values 
as we could for the Hubbard model. In fact, 6 is larger than 4 and continues to increase 
for finite densities and increasing negative interaction strength s ^ — 1^. A plot of as a 
function of the density do for different values of the interaction strength s is given in Fig. |^. 
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For bose statistics, Dq and D_i are restricted to integer values. The correlators of 
diagonal operators, i.e. the density-density correlator C„ and the longitudinal spin-spin 
correlator are independent of statistics, and so only the correlators C^, and C, 



sing 



change. We find for their asymptotics 



1 

\x — ivotJ \x — iv-it J 



+ 



Re 



\x + tV-it\^^^^> 



\x 



-ifot|''/4|x + if„it|2+^'/(i+^) 



Re 



\x — iv^it 



\x + ivot\^/'^+'^/^\x + iv-it 



1/(1+^) 



Re 



Aie 



\x — ivot 



(53) 
(54) 



(55) 



(56) 



V. ASYMPTOTICS OF THE CORRELATION FUNCTIONS FOR THE BOUND 

CASE 



Due to the restriction ([20| ) on the quantum numbers of a given state, Di is an integer 
or half-odd integer depending on the parity of AA^^i for both bose and fermi statistics of the 
particles, i.e. 



(mod 1) 



(57) 



This selection rule is just the same as the case of one-component bosons, and so we find for 
the asymptotics of the pair density correlator 

- {vitf 



Cp{x, t) -dj^ A 



+ (t;it)2]2 

+A2 cos{2kfX + (pi) 



\x + ivit\^ 



(58) 
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and for the pair field correlator 
C^{x,t) ~ Ai- 



^ rRe 



X + ivit 



(59) 



\x + ivit\<^+^/^ 

Here we again defined an exponent 6 = 2^^. Following equation (0), we can calculate 
^1 from the Fermi velocity of pairs fi. In Fig. ^ we plot the lines of constant in the 
[do, s) plane. Note that the value of ^(^i) at zero density is given by 8(2), whereas for finite 
densities and vanishing interaction strength s — 0+, we have ^ ^ 4 (^i ^ V^)- A plot of 
^1 as a function of the density do for different values of the interaction strength s is given in 
Fig. |. 

VI. THE NON-INTERACTING TWO-COMPONENT SYSTEM 

At s = 0, the system reduces to a non-interacting two-component gas and we may ex- 
pect a certain continuity in the behavior of the correlators at this point. Indeed, as s — 0~, 
the two Fermi velocities vq and f_i both approach the Fermi velocity of a non- interacting 
single-component model, i.e. vo{s 0^) = f_i(s 0^) = 7rc/o/2. Consequently, the cor- 
relation functions of the bosonic (fermionic) system reduce to the correlation functions of a 
non-interacting bose (fermi) system with two components, i.e. with half the one-component 
fermi momentum. Using the language of conformal field theory, we can thus describe the ex- 
citations of the non-interacting two-component gas by a c = 2 generalized Gaussian model.§ 

From the expression of the dressed charges and ^i, we see that = \^o^- As s ^ 0+, 
the Fermi velocity of the pairs goes to the Fermi velocity of a one-component free bose gas 
with doubled particle mass, i.e. vi{s 0^) = 'Kdi/2 = ^vq{s 0^). Therefore, we expect 
6i = 6q at s = and this is indeed true as shown above. Furthermore, the free energy of 
the system should be uniquely specified at s = 0. Following (^3|) we may write the finite 
temperature corrections for the unbound 

F T ~ F T = ---— + - , 60 

6 \v^i Vol 



whereas for the bound case we have 

F{T) ^ F{T = 0) - (61) 

As predicted, these two equations are in agreement at s = and identical to the free energy 
of a non-interacting c = 2 system. 

The bound pairs for s > are singlets. Therefore we might expect that the pair field 
correlator (^6|) becomes identical to the singlet pair correlators ( ^4|) and ( ^6l) of the unbound 



s — > 0. However, \E'^ creates pairs with characteristic length scale 1/s and not just 
two particle wave functions. Thus, the pair wave functions include a normalization factor 
^/s. As s ^ 0"*", the leading terms of the conformal expansion (|59|) consequently vanish 
and higher order terms become important. It should therefore come as no surprise that the 
expansions (H), (0) and (|59D do not agree at s = 0. 
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FIGURES 

FIG. 1. Lines of constant universal behavior for the unbound case. Contours of constant value 

of the dressed charge hi the {do, s) plane arc shown. The lines represent increments of .2 starting 
from ^0 = 1-0 at do = up to = 1-8- The dashed line correspond to the value = of a 
non-interacting system. 

FIG. 2. Plot of 6 as function of particle density do for various values of interaction strength s 
for the unbound case. 

FIG. 3. Lines of constant universal behavior for the bound case. Contours of constant value of 
the dressed charge in the {do,s) plane are shown. The lines represent increments of .2 starting 
from ^1 = 2.0 at do = down to = 1.2. 

FIG. 4. Plot of 9 as function of particle density do for various values of interaction strength s 
for the bound case. 
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